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ABSTRACT 


Many  applications  lead  to  nonlinear,  parameter  dependent  equations 
H(y,t)  *  yQ,  where  H:  Y  *  T  Y,  yQ  e  rge  H,  and  the  state  space  Y  is 
infinite-dimensional  while  the  parameter  space  T  has  finite  dimension. 

The  case  dim  T  =  1  is  of  special  interest  in  connection  with  continuation 
methods.  For  this  case,  a  general  theory  is  developed  which  provides  for  the 
existence  of  solution  paths  of  a  rather  general  class  of  Such  equations  and 
of  their  finite-dimensional  approximations,  and  which  allows  for  an  assessment 
of  the  error  between  these  paths.  A  principal  tool  in  this  analysis  is  the 
theory  of  nonlinear  Fredholm  operators.  The  results  cover  a  more  general 
class  of  operators  than  the  mildly  nonlinear  mappings  to  which  other  approaches 
appear  to  be  restricted. 


On  the  Discretization  Error  of  Parametrized 


by 

James  P.  Fink^  and  Werner  C.  Rheinboldt^ 


1 .  Introduction 

Nonlinear  parameter-dependent  equations  of  the  form 

H(y,t)»y0  (1.1) 

arise  in  connection  with  many  equilibrium  problems  in  science  and  engineering. 

Here  y  and  t  vary  in  some  state  space  Y  and  parameter  space  T,  respectively, 
H  is  a  mapping  from  Y  x  1  into  Y,  and  yQ  is  a  given  point  in  the  range  of 
H.  Under  appropriate  conditions,  the  set  of  all  solutions  (y.t)  of  (1.1)  forms 
a  manifold  in  Y  x  T  and  we  wish  to  analyze  the  properties  of  this  manifold. 

In  most  cases,  the  state  space  Y  is  infinite-dimensional  and  the  parameter 
space  T  finite-dimensional,  although  infinite-dimensional  parameter  spaces  have 
been  considered  as  well  (e.g.,  see  [1]).  Thus,  for  practical  computations,  there 
is  a  need  for  introducing  finite-dimensional  approximations  of  the  equation  (1.1). 
This  leads  to  questions  about  the  relation  between  the  properties  of  the  solution 
manifolds  of  the  original  and  the  discretized  equations  and  the  errors  introduced 
by  the  approximation.  Relatively  little  has  been  done  so  far  in  this  area.  With¬ 
out  any  claim  of  completeness,  we  mention  here  only  the  articles  [2],  [5],  [7], 
[12],  [15]  which,  in  different  settings,  address  the  a  priori  estimation  of  the 

1)  This  work  was  supported  in  part  by  the  U.S.  Air  Force  Office  of  Scientific 
Research  under  Grant  80-0176,  the  National  Science  Foundation  under  Grant 
MCS-78-05299  and  the  Office  of  Naval  Research  under  Contract  N-0014-80-C-0455. 

2)  Department  of  Mathematics  and  Statistics,  University  of  Pittsburgh,  Pittsburgh, 
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approximation  errors  in  the  case  of  a  one- dimensional  parameter  space.  In 
[3],  some  a  posteriori  estimates  of  these  errors  are  considered  for  certain 
two-point  boundary-value  problems. 

The  restriction  to  one-dimensional  parameter  spaces  is  not  as  severe  as 
it  might  appear.  In  fact,  all  available  computational  methods  for  analyzing 
solution  manifolds  of  (finite-dimensional)  parametrized  equations  consist 
principally  of  some  form  of  continuation  process  for  the  trace  of  paths  on  the 
manifold.  Such  paths  are  the  solution  manifolds  of  some  equation  involving  a 
one-dimensional  parameter. 

In  this  paper,  we  present  a  general  theory  which  provides  for  the  existence 
of  solution  paths  of  a  rather  general  class  of  nonlinear  equations  with  a  one¬ 
dimensional  parameter,  as  well  as  of  their  finite-dimensional  approximations, 
and  which  allows  for  an  assessment  of  the  error  between  these  paths.  A  principal 
tool  in  this  analysis  is  the  theory  of  nonlinear  Fredholm  operators.  As  in  [5], 
[7],  the  error  estimates  are  derived  from  a  form  of  the  implicit  function  theorem. 
But  our  results  cover  a  more  general  class  of  operators  than  the  mildly  nonlinear 
mappings  to  which  the  approaches  in  [5],  [7]  are  restricted. 

2.  Solution  Manifolds  and  Local  Parametrizations 

In  order  to  avoid  repetition,  we  shall  assume  throughout  this  paper  that 
the  following  information  has  been  given: 

(i)  X,Y  real  Banach  spaces; 

(A) 

(ii)  F:  Ecx  +  Y  a  mapping  of  class  Cr ( E ) ,  r  >_  1,  on  some 
open  subset  E  e  X. 

In  finite  dimensions,  an  equation  involving  an  m-dimensional  parameter  corresponds 
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to  the  case  dim  X  -  dim  Y  =  m.  The  analog  in  infinite  dimensions  is  the 
requirement  that  ker  DF(x)  has  dimension  m  for  the  points  x  under 
consideration.  More  specifically,  we  are  interested  in  regular  points  of 
F,  that  is,  those  x  e  E  for  which  DF(x)  maps  X  onto  Y.  Hence,  we 
assume  that  the  m-regularity  set 

Rm(F)  =  {x  e  E:  dim  ker  DF{x)  =  m,  DF(x)X  =  Y}  (2.1) 

is  nonempty  for  some  m  >_  1.  Then  the  restriction  of  F  to  R  (F)  is  a  non¬ 
linear  Fredholm  operator  of  index  m  (e.g.,  see  [4]),  and  this  observation  is 
at  the  center  of  our  analysis. 

Lenina  2.1:  The  sets  R  (F),  m  >  1 ,  are  either  empty  or  open  in  X. 

Proof:  Suppose  that  x  e  Rm ( F ) .  Then  by  the  continuity  of  DF  and  a  basic 
perturbation  property  of  Fredholm  operators  (e.g.,  see  [13,  p.  115])  there  exists  a 
6  >  0  such  that  DF(x+£)  is  again  a  Fredholm  map  of  index  m  for  all  E,  e  X 
with  1 15|  |  <  <5.  Moreover,  we  have  dim  ker  DF(x+£)  £  dim  ker  DF(x),  whence 
dim  coker  DF(x+C)  =  0  and  dim  ker  DF(x+£)  =  m. 

We  are  interested  in  examining  the  solution  set  of  the  equation 

F(x)  =  yQ,  x  e  Rm(F),  (2.2) 

for  a  fixed  yQ  e  F(Rm(F)),  that  is,  the  set 

W  m  »2-3> 


For  any  xQ  c  Rm(F),  there  exist  closed  subspaces  Vc  X  such  that 
X  =  V  ©  leer  DF(xq).  For  any  such  choice  of  V,  the  restriction  DF(xQ)|V 
is  an  isomorphism  between  V  and  Y,  and  its  inverse 

Ay  -  (DF{x0)|V)_1  e  L(Y,V)  (2.4 

is  a  right  inverse  of  DF(xQ).  With  this,  equation  (2.2)  may  be  written  in 
the  form 

F(t  +  Avy)  =  yQ,  t  e  ker  DF(xq),  ytY,  t  +  Ayy  e  RJF),  (2.5 

which  corresponds  to  (1.1)  with  Y  as  the  state  space  and  T  =  ker  DF(xq) 
as  the  m-dimensional  parameter  space. 

The  basic  result  about  the  solution  set  (2.3)  may  now  be  phrased  as 
follows: 

Theorem  2.2:  Suppose  that  the  information  (A)  is  given  and  that  R  (F)  is  non¬ 
empty  for  some  m  >  1.  Then,  for  any  yQ  e  F(Rm(F)),  the  regular  solution  set 
Mm(y0)  is  a  nonempty,  relatively  open,  m-dimensional  Cr-manifold  in  X. 

Proof:  The  proof  proceeds  along  the  lines  of  the  finite-dimensional  analog 
given  in  [10],  Let  xQ  e  Mm(yQ)  and  V  be  a  subspace  of  X  such  that 
X  *  V  ©  ker  DF(xq).  With  the  corresponding  projection  P:  X  -*•  ker  DF(xq),  we 
define  the  mapping  G:  E  -*■  Y  *  ker  DF(xq)  by  G(x)  =  (F(x),Px)  for  x  e  E. 

Then  DG(x)£  =  (DF(x)£,P£).  (  e  X,  and  DG(xo)  is  nonsingular.  Hence,  the 
inverse  function  theorem  applies  and  G  maps  some  neighborhood  ilex  of  xQ 
Cr-diffeomorphically  onto  a  neighborhood  S  of  (y0,PxQ).  Evidently  now,  G 
maps  Mm(yQ)  n  U  Cr-diffeomorphically  onto  ({yQ} *  ker  DF(xQ) )  n  7  and,  since 
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{yQ}  x  ker  DF(xQ)  Is  Cr-diffeomorphic  to  IRm,  the  result  follows. 


As  mentioned  in  the  introduction,  the  computational  procedures  for  analyzing 
the  solution  manifold  Mm(yQ)  of  (2.2)  consist  principally  of  methods  for 
approximating  paths  on  this  manifold.  Such  a  path  is  defined  as  the  solution 
manifold  of  some  reduced  equation  defined  by  a  Fredholm  operator  of  index  1 
on  its  regularity  set.  Therefore,  from  now  on  we  restrict  ourselves  to  the 
case  when  (A)  is  given  and  R-j(F)  is  nonempty.  For  ease  of  notation,  we  will 
write  R(F)=R](F)  and  M(yQ)  =  M] (yQ) . 

We  consider  first  the  question  of  the  choice  of  suitable  local  parametrizations 
of  the  one-dimensional  manifold  M{yQ).  It  turns  out  that  a  possible  choice 
corresponds  to  a  typical  approach  used  in  many  continuation  procedures. 

In  line  with  the  earlier  discussion  leading  to  (2.5),  a  local  parametrization 
of  M(yQ)  at  a  given  point  xQ  e  M(yQ)  is  defined  as  a  triple  {V,A,zq}  con¬ 
sisting  of  a  closed  subspace  V  of  X,  a  linear  operator  A  e  L(Y,V),  and  a  point 
zQ  e  X  such  that 


(i )  V  n  ker  DF(xq)  =  {0} ; 

(ii)  A  is  an  isomorphism  of  Y  onto  V; 

(iii)  zc  i  V; 

( i v)  X  =  V  ©  T0,  TQ  =  span  {zq}. 


(2.6) 


As  indicated  in  Figure  1,  we  consider  the  family  of  parallel  linear  mani¬ 


folds 
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xo  +  tz  o  +  V*  t£lRl*  (2.7) 

The  conditions  (2.6)  ensure  that  these  manifolds  are  transverse  to  ker  DF(xQ) 
and  hence  that,  locally  near  t  =  0,  the  solution  manifold  may  be  parametrized 
in  terms  of  t.  This  is  the  content  of  the  following  result: 

Theorem  2.3:  Let  (A)  be  given,  R(F)  f  <t>,  and  {V,A,zq}  a  local  para- 

metr-'zation  of  M(yQ)  at  a  point  xQ  e  M(yQ).  Then  there  exists  an  open 
interval  Jc]r',  0  e  J,  an  open  neighborhood  U  e.x  of  xQ,  and  a  unique 
Cr-map  n-‘  J  Y,  such  that 

M(yo )  H  U  =  {x  e  X:  x  =  xQ  +  tZQ  +  An(t),  t  e  J).  (2.8) 

Proof:  Let  B:  1R1  *  Y  -*  X  denote  the  affine  mapping  B(t,y)  =  xQ  +  tzQ  +  Ay, 
t  e  IR1 ,  y  e  Y,  and  define  H:  8^"^(E)  +  Y  by 

H(t,y)  =  F(B(t,y) )  - y0,  (t.y)  e  B(_1)(E).  (2.9) 

Then  H  is  of  class  Cr  on  B^”^(E)  and  H(0,0)  =  0.  Moreover,  since  A 
is  an  isomorphism  of  Y  onto  V  and  DF( xQ )  maps  V  one-to-one  onto  Y,  it 
follows  that  DyH(0,0)  =  DF(xq)A  is  an  isomorphism  of  Y.  Therefore,  by  the 

implicit  function  theorem,  there  exists  an  open  interval  J  c  IR^ ,  0  e  J,  a 

neighborhood  S  «  Y  of  y  =  0,  and  a  unique  Cr-function  q:  J  -*■  S  such  that 

H(t,n( t) )  =  F(xq  +  tZQ  +  An(t))  -  yQ  =  0,  t  e  J. 

With  U  =  B(d  x  S),  the  uniqueness  of  n  ensures  that  (2.8)  holds. 
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In  the  standard  continuation  methods,  zQ  defines  the  predictor  line  and 
the  corrector  produces  the  step  An(t)  from  the  predicted  point  xQ  +  tzQ 
to  the  point  xQ  +  tzQ  +  An(t)  on  M(yQ).  A  particular  local  parametrization 
is  obtained  if  we  choose  any  nonzero  vector  zQ  e  ker  DF(xQ )  and  the  map  Ay 
of  (2.4)  for  any  closed  subspace  V  c  X  such  that  X  =  V  ©  ker  DF(xQ).  This 
choice  corresponds  essentially  to  the  pseudo-arclength  parametrizations  proposed 
by  various  authors  (e.g.,  see  [8]). 

The  question  arises  as  to  how  far  a  particular  local  parametrization  {V,A,zq} 
can  be  extended.  For  this,  note  that  the  set 

A  =  {x  e  E:  DF(x)A  is  an  isomorphism  of  Y}  (2.10) 

is  certainly  nonempty  and  open.  With  this  observation,  a  generalization  of  a 
result  in  [9]  may  be  phrased  as  follows: 


Theorem  2.4:  Under  the  conditions  of  Theorem  2.3,  let  M  c  X  denote  the 
maximal  connected  subset  of  M(yQ)  rt  A  which  contains  xQ.  Then  there  exists 

•  1  r 

an  open  interval  J  e  ]R  ,  0  e  J  ,  and  a  C  -function  n  :  J  Y  such  that 

■  OO  00  00  oo 


M  =  {x  e  X:  x  =  xn  +  tzn  +  An  (t),  t  e  J  }. 

oo  0  0  00  00 


Proof:  Since  X  =  V  ©  T  ,  each  x  e  X  can  be  written  uniquely  in  the  form 
x  =  xQ  +  tzQ  +  Ay  for  (t,y)  e  *  Y.  Thus,  each  x-|  e  may  be  expressed 

uniquely  as  x-j  =  xQ  +  t-j zq  +  Ay-j  for  (t^,y^)  e  B^"^(E),  where  B  is  defined 
in  the  proof  of  Theorem  2.3.  Let  H:  B^’^(E)  +  Y  be  given  by  (2.9).  Then, 
for  any  x1  e  M^,  the  derivative  DyH(t1,y1)  =  DFCx^A  is  an  isomorphism  of  Y, 
and  we  may  apply  the  implicit  function  theorem  to  ensure  the  existence  of  an 
open  interval  J  c  1R1 ,  t,  e  J,  ,  a  neighborhood  Sv  c  Y  of  y, ,  and  a  unique 

X-j  I  X-j  X^  I 
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C  -function  n  :  J 
X1  X1 


such  that  H(t,n  ( t ) )  =  0  for  t  e  J  .  Con- 


'1 


1 


A1 

sider  the  projection  it:  1R*  *  Y  IR1  and  its  restriction  to 

Then,  for  each  e  M^,  and  hence  each  (t^  ,y1 )  e  N^,  there  exists  an  open 

interval  Jx  e  IR^  containing  t-j  such  that  tt^  =  irjW^  is  a  one-to-one 

mapping  of  )  onto  J  .  In  fact,  we  have 


=  (t,n  ( t) ) ,  t  e  J  . 

X1  X1 

r  1 

Thus,  nm  is  a  local  C  -diffeomorphism  of  to  1R  .  A  classical  result  of 
topology  states  that,  as  a  local  homeomorphism  from  the  connected  Hausdorff 
space  W  to  !r\  h  must  be  a  homeomorphism  from  N  onto  J  =  tt  (W  ). 
Now  with  the  other  projection  4>:  k'  *  Y  Y  we  may  define  nro:  -*■  Y  by 

n  =  0  °  if1  and  the  result  follows. 

co  T  oo 


The  condition  that,  for  x  e  MQ,  the  mapping  DF(x)A  is  an  isomorphism 
of  Y  means  geometrical ly  that  a  local  parametrization  {V,A,zq}  is  valid  for 
the  segment  of  the  solution  path  M(yQ)  between  the  points  x_  and  x+  closest 
to  xQ  at  which  x_  +  V  and  x+  +  V  are  tangent  to  M(yQ)  (see  figure  1).  In 
standard  terminology,  these  are  the  closest  limit  points  of  the  path  with  respect 
to  the  decomposition  X  =  V  ©  span(zo). 


3.  Finite-Dimensional  Approximations 

In  this  section,  we  turn  to  the  formulation  of  suitable  approximate  problems 
for  (2.2).  In  many  applications,  one  has  X  =  Y  *  1R1 ;  that  is,  a  particular 
component  of  X  is  identified  as  a  parameter  and  only  the  complementary  component 
Y  needs  to  be  discretized.  We  generalize  this  situation  by  assuming  that  a 
splitting 

! 
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X  =  Z  ©  T,  dim  T  *  1 


has  been  given  together  with  an  operator  which  relates  Z  to  Y:  namely. 


(3.1} 


Q  e  L(X,Y) ,  T  =  ker  Q, 


Q | Z  e  L(Z,Y)  is  an  isomorphism. 


(3.2) 


As  in  the  previous  section,  suppose  that  our  basic  data  (A)  are  given  and 
that  R( F)  =  R-|(F)  is  nonempty.  A  discretization  of  the  problem  is  specified 
by  a  set  {P^}  of  projections  e  L(Y)  of  finite  rank,  indexed  by  positive 
h  >  0,  such  that 


lim  P.y  =  y,  ye  Y. 
h-0  n 


(3.3) 


With  the  families  of  subspaces 


Yh =  V-  V  wfX,  xh  =  zh©T, 


h*  h  h 


(3.4) 


the  approximate  problem  is  then  given  by 


where 


Fh(x)  =  yoh,  x  e  R(Fh).  yoh  -  Phy0, 


Fh:  Eh  c  xh  "*■  Yh’  Fh^x^  =  PhF^’  x  E  =  E  n  Xh‘ 


(3.5) 


(3.6) 


For  ease  of  discussion,  we  call  the  information  given  by  the  projections  {P^}, 
the  subspaces  (3.4),  and  the  approximate  equations  (3.5)  and  operators  (3.6)  a 
basic  approximation  of  the  problem  (2.2). 


't  *.j 


4- 


The  general  question  we  address  in  this  and  the  next  section  is  whether 
the  solution  manifolds  of  the  discretized  equations  (3.5)  approximate  the 
solution  manifold  M(yQ)  of  the  original  equation  (2.2)  when  h  tends  to 
zero.  More  specifically,  in  this  section  we  consider  the  existence  of 
solutions  of  (3.5)  and  take  up  the  matter  of  error  estimates  in  the  subsequent 
section. 

For  the  analysis  of  the  approximate  problems,  it  is  convenient  to  extend 
the  discrete  operators  (3.6)  to  all  of  E  e  x  as  follows: 

/. 

Fh:  E  -  Y,  Fh(x)  =  (I-Ph)Qx  +  Ph(F(x)-yQ),  x  e  E,  (3.7) 

1 

■ ! 

where  I  denotes  the  identity  on  Y.  Clearly,  Fh  is  a  Cr-mapping  and  the  \ 

following  properties  hold:  i 

Proposition  3.1: 

A 

(i)  F^(x)  =  0  ^or  x  e  E  ^  and  °Hly  ^  x  e  Eh  and  F^(x)  =  yQh; 

(1i )  °Fh(x)xh  cYh,  xeE; 

(iii)  ker  DFh(x)c  X^,  x  e  E; 

(iv)  DFh(x)  e  L(X,Y)  is  a  Fredholm  operator  of  index  1  for  all  x  e  E;  1 

(v)  P^DF(x)Xh  =  Yh  for  some  x  e  E  implies  that  x  e  R(Fh). 

Proof:  The  proofs  of  (i),  ( i i ) ,  and  (iii)  are  straightforward.  For  example,  if 
F^(x)  =  0  for  some  x  e  E,  then  Qx  =  PhQx  -  P^( F(x)-yQ )  e  Y^  implies  that 
x  e  Xh  and  hence  that  Fh(x)  =  PhF(x)  =  Fh(x)  -  (I-Ph)Qx  +  PhyQ  =  yQh.  The 
converse  of  (i),  as  well  as  (ii)  and  (iii),  follow  analogously.  The  property  j 

. . . .  . 
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(iv)  is  a  consequence  of  the  fact  that  DF^x)  =  Q  -  PhQ  +  PhDF(x)  is  a  compact 

perturbation  of  the  Fredholm  operator  Q  e  L(X,Y)  of  index  1.  Since 

DF^ (x ) | Xh  =  PhDF(x)Xh,  we  obtain  (v)  from  (iii)  and  (iv)  if  only  it  can  be 

A 

shown  that  dim  ker  DFh(x)  =  1.  But,  since  by  assumption  rge  Ph DF( x ) ) X^  =  Yh 
and  dim  Xh  =  dim  Y^  +  1 ,  this  follows  immediately  from  dim  ker  PhDF(x)|Xh  + 
dim  rge  P^DF(x)|X^  =  dim  X^. 

Clearly,  any  comparison  of  the  solution  manifolds  of  (2.2)  and  (3.5)  must  be 
done  locally.  Let  xq  e  M(yQ)  be  given  and  suppose  that  (V,A ,zQ}  is  a  local 
parametrization  of  M(yQ)  at  xQ.  This  parametrization  must  relate  in  a  suitable 
way  to  the  basic  discretization  introduced  above,  and  the  compatibility  condition 
may  be  expressed  in  a  variety  of  forms.  Rather  than  go  into  a  detailed  discussion 
of  various  possible  equivalent  definitions,  we  introduce  here  simply  the  following 
technical  condition  which  is  relatively  easy  to  verify  in  many  practical  situations 

(C)  | J DFh ( xq )Ay | |  ^  Y | |y | | ,  y  e  Y,  for  all  sufficiently  small  h  >  0. 

Here  y  >  0  is  a  constant  independent  of  y  and  h. 

Compatibility  conditions  of  this  form  typically  arise  in  the  consideration  of 
projection  methods.  For  example,  condition  (C)  is  related  to  the  notion  of  stable 
convergence  of  linear  operators  defined  and  discussed  in  [14]. 

The  condition  (C)  can  be  enforced  in  a  number  of  ways.  For  example,  by  a 
judicious,  although  somewhat  restrictive,  choice  of  the  discretization  and  pa¬ 
rametrization,  DFh(xQ)A  turns  out  to  be  the  identity  on  Y  and  (C)  holds  with 
y  =  1  for  all  h  >  0.  As  discussed  in  Section  7,  this  case  includes  the  so- 
called  reduced-basis  technique  (e.g.,  see  [11]).  On  the  other  hand,  as  noted 
at  the  beginning  of  this  section,  in  many  applications  the  formulation  of  the 


problem  includes  a  natural  splitting  of  the  space  X,  and  then  (C)  becomes 
in  essence  a  condition  on  F  if  we  are  not  prepared  to  restrict  the  dis¬ 
cretizations  and  parametrizations.  In  particular,  we  show  in  Sections  5  and 
6  several  classes  of  operators  F  for  which  (C)  is  easily  verified. 

As  a  first  consequence  of  (C),  we  obtain  the  following  result: 

Proposition  3.2:  If  (C)  holds,  then  PhDp(*0)Xh  =  \  and  x0  e  R(Fh)  for 
all  sufficiently  small  h  >  0. 

A 

Proof:  Evidently,  (C)  implies  that  DF^(xq)A  e  L(Y)  is  an  isomorphism.  Hence, 
for  y^  e  Y^,  there  exists  a  y2  £  V  such  that  DFh(xo)Ay2  =  y-j ,  and  therefore 

QAy2  =  PhQAy2  “  phDF^xo^Ay2  +  yl  e  Yh‘  From  this  follows  that  Ay2  e  xh’ 
(I-ph)QAy2  =  °»  and  PhDF^x0^2  =  yl  ‘  The  second  assertion  is  now  a  direct 
consequence  of  Proposition  3.1  (v). 

In  order  to  prove  the  existence  of  solutions  of  the  approximate  problems 
(3.5),  we  make  use  of  a  generalized  form  of  the  inverse  function  theorem.  As 
noted  in  [5],  the  usual  proof  of  this  theorem  provides  a  Lipschitz  condition 
for  the  inverse  function.  In  our  terminology,  the  result  mav  be  formulated  as 
fol 1 ows : 

Theorem  3.3:  Let  (A)  be  given,  where,  at  some  xo  e  E,  (i)  DF(xQ)  e  L(X,Y)  is 
an  isomorphism  onto  Y  with  1 1 0F(x  ) |  }  <  1/y,  and  (ii)  there  exists  a  6  >  0 
such  that 

|  | DF(x)-DF{xq)  |  |  <  J,  x  e  8{xo,5)  =  {x  e  <:  |  |x-xQ|  |  <  S)CE. 

Then  there  exists  a  unique  Cr-map  G:  8(F(xq)  ,y<$/2)  8(xQ,6)  for  which 
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F(G(y) )  =  y  for  all  y  e  S(F(xq) ,y«/2)  and 

I  |G(y1  )-G(y2)  1 1  <  ~  |  |y-, -y2 1 1  •  y^  e  B(F(x0),y6/2). 

In  order  to  prepare  the  way  for  the  application  of  this  theorem,  let 
B:  Y  -*■  X  be  the  affine  mapping  given  by  By  =  xQ  +  Ay,  y  e  Y,  and  define 
the  operators 

Hh:  B(_1,(R(F))  +  Y,  Hh(y)  =  Fh(By),  y  e  B*"1 }(R(F) ) .  (3.8) 

Clearly,  is  of  class  Cr  and  DH^Cy )  =  (I-Ph)QA  +  PhDF(xQ+Ay)A. 

We  wish  to  apply  Theorem  3.3  to  The  two  conditions  of  the  theorem 

are  verified  in  the  following  lemma: 

Lemma  3.4:  Let  (C)  hold.  Then  DHh (0)  e  L(Y)  is  an  isomorphism  of  Y  with 
| | DH^(O) | |  £  1/y  for  all  sufficiently  small  h  >  0.  Moreover,  there  exists 
a  6  >  0,  independent  of  h,  such  that 

1 ) DHh(y)-DHh(0) | |  <  J  whenever  ||y||<6.  (3.9) 

Proof:  The  first  part  is  an  immediate  consequence  of  (C).  In  order  to  prove 
(3.9),  note  first  that,  from  (3.3)  and  the  uniform  boundedness  principle,  it 
follows  that 

|  (Ph(|  a  for  all  sufficiently  small  h  >  0  (3.10) 


with  a  constant  a  >  0  which  does  not  depend  on  h.  Hence,  we  find  that 
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||DHh(y)-DHh{0)||  =  ||Ph(DF(x0+Ay)  -  DF(x0))A|| 


<  a|| A||  ||DF(x0+Ay)  -  DF(x0)||, 
and  (3.9)  is  a  direct  consequence  of  the  continuity  of  DF. 


With  this,  we  may  now  apply  Theorem  3.3  to  Hh  for  sufficiently  small 
h  >  0.  This  ensures  the  existence  of  a  unique  Cr-function 


Gh:  8o  =  8<(I-Ph>QV*/2>  "  B{0'6) 

for  which 

(i)  Hh(Gh(y))  =  y,  y  £  Bq, 

(ii)  ||Sh(y1)-Gh(y2)||  <  |  ||y1-y2||,  y, ,y2  e  BQ. 


(3.11) 


(3.12) 


Since  y  and  6  are  independent  of  h  and  (I-P^QXq  0  as  h -*■  0,  we  see 

that  0  e  Bq  for  all  sufficiently  small  h  >  0.  Hence,  yh  =  Gh(0)  satisfies 

Hh(yh)  =  0,  and  from  (3.12)  (ii)  with  y^  =  0  and  y2  =  (I-Ph)QxQ  it  follows 
that  | |yh| |  <  (2/y )  | | (I-Ph)QxQ| | .  Hence,  for  xQh  =  xQ  +  Ayh,  we  obtain 
Fh(xQh)  =  Hh^yh^  =  0  as  well  as  x0h  x0  when  h  "*■  0*  Moreover,  Proposition 

(3.1)  (i)  implies  that  xQh  £  Eh  and  F^(xQfl)  =  yoh‘  *n  °ther  words,  xQh  is 

the  desired  solution  of  the  approximate  problem  (3.5)  corresponding  to  xQ.  These 
solutions  have  several  interesting  properties: 


Proposition  3.5:  If  (C)  holds,  then  xQ^  £  R( F^ )  for  all  sufficiently  small  h. 
Proof:  Since  R(F)  is  open,  we  certainly  have  xQh  e  R(F)  for  all  sufficiently 

A 

small  h.  Moreover,  by  Proposition  3.2,  we  also  know  that  xQ  e  R(Fh).  But,  in 
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order  to  conclude  that  the  x^  ultimately  belong  to  R(Fh),  we  need  to  use 
the  uniformity  implied  by  (C).  For  this,  we  write 

Dph(xoh)A  *  Dph(xo)A  +  <0fh<xoh>-Dph<xo»A 

and  note  that,  for  all  sufficiently  small  h  >  0,  it  follows  from  Lemma  3.4 
that  DFh(xo)A  e  L(Y)  is  an  isomorphism  with  1 1  ( DF^(xQ)  A) 1 1  <_\/y  and, 
moreover,  that,  by  (3.9)  with  y  =  y^. 


li(Dfh<xoh>-rfh(xo>>All  i  }■ 

Therefore,  °^^xoh^A  e  L(Y)  is  also  an  isomorphism  of  Y  and  thus  D^(xoh) 

A 

must  map  X  onto  Y.  Since,  by  Proposition  (3.1)  (iv),  is  a  Fredholm 

operator  of  index  1,  the  result  follows. 

Proposition  3.6;  If  (C)  holds,  then  V  n  ker  DF^(xQh)  =  {0}. 

Proof;  In  the  proof  of  Proposition  3.5,  we  saw  that  DFh(xQh)A  e  L(Y)  is  an 
isomorphism  of  Y  for  sufficiently  small  h.  Hence,  if  D^xoh^v  =  0  for 

A 

some  v  e  V,  then  v  =  Ay  for  some  y  e  Y,  and  therefore  DFh(xQh)Ay  =  0, 
whence  y  =  0  and  thus  also  v  =  0. 

We  collect  our  results  up  to  this  point  in  the  following  theorem: 

Theorem  3.7:  Let  (A)  be  given  and  R(F)  t  <(>.  Suppose  that  a  basic  approximation 
has  been  chosen  and  that  {V,A,zQ}  is  a  local  parametrization  of  M(y0)  at 
x0  c  M(yQ),  where  the  compatibility  condition  (C)  is  satisfied.  Then,  for  all 
sufficiently  small  h  >  0,  the  approximate  problems  (3.5)  have  solutions 


V-  Mfr' 
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xoh  e  R(Fh}  for  wh1ch  "J"  xoh 

h-KJ 


>-i 


of  the  form  { Vh , ( PhDF( xQh ) | Vh) 
of  (3.5)  at  xQ^,  Theorems  2.2, 


xQ.  Moreover,  for  any  local  parametrization 

,zoh}  of  the  solution  manifold  F^-1^(yoh)  n  R( ) 
2.3  and  2.4  apply. 


4.  Error  Estimates 

Theorem  3.7  ensures  the  existence  of  a  solution  segment  for  the  local 
approximate  problem  (3.5)  under  the  same  general  conditions  on  the  operator  F 
needed  to  establish  the  existence  of  a  solution  segment  for  the  full  problem 
(2.2).  If  error  estimates  are  desired,  however,  then  additional  smoothness 
conditions  are  required  for  F.  Moreover,  the  parametrizations  of  the  two 
curves  cannot  be  chosen  independently,  but  must  be  such  as  to  allow  for  a  specific 
association  of  comparable  points. 

For  the  derivation  of  the  error  estimates,  we  use  a  form  of  the  implicit 
function  theorem  formulated  in  [5]  which  follows  from  Theorem  3.3. 

Theorem  4.1 :  Let  X,Y,Z  be  Banach  spaces,  and  $:  S  «=  X  -*■  Y  a  Lipschitz- 
continuous  mapping  with  Lipschitz  constant  >  0  on  the  subset  S  of  X. 

Let  F:  EeXxY  +  Z  be  a  Cr-mapping  on  a  neighborhood  E  of  S  x  3>(S)  such 
that  (i)  for  each  x  e  S,  D  F(x,4>(x))  is  an  isomorphism  of  Y  onto  Z  for 
which  sup{|  |DyF(x,$(x))_1  |  | :  x  e  S}  _<  yq  and  sup{  |  |  DxF(x,$(x) )  |  | :  x  e  S}  <_  y1  , 
and  (ii)  there  exists  a  monotonically  increasing  function  l:  rcj  -*■  ]R^  for 

|  |DF(x,y)-DF(x0,$(x0))  ]  |  £  £(p)  I  !  ( x,y ) - ( xQ, «^( xQ) )  |  j ,  xq  e  S, 

(x.y)  e  8((x0,*(x0)),p). 


which 
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Then  there  are  constants  a  >  0  and  b  >  0,  depending  only  on  and 

£,  such  that,  for  sup{  |  |F(xo,$U0)  1 1 :  xQ  e  S}  <  e  =  a/2,  there  exists  a 

Cr-function  G:  8(x.,a)  Y  which  satisfies 

x  eS  0 
0 

F(x,G(x) )  =0,  x  c  U  8(x  ,a), 
xeS  0 

°  (4.1) 

G(8(x0,a))  c  8(*U0),b),  xQ  e  S, 

and 

||S(x)-*(x0)||  <  C0(||x-x0||  +  ||F(x0,*(x0))||),  xQ  e  S,  x  e  8(xo>a), 

(4.2) 

where  the  constant  CQ  >  0  depends  only  on  y0  and  y, . 

In  order  to  set  up  a  situation  in  which  we  can  use  this  theorem,  we  continue 
to  assume  the  conditions  of  Theorem  3.7.  By  Theorem  2.3,  the  original  problem 
(2.2)  has  a  solution  segment  x(t)  =  xQ  +  tzQ  +  Ari(t)  defined  for  t  in  an 
open  interval  jclR1  containing  t  =  0.  By  the  continuity  of  OF,  there  exists 
a  nontrivial  compact  subinterval  JQ  c  «i,  0  e  JQ,  such  that 

llDF(xo  +  t2o  +  Mt))  -  DF(x0)||  £  8  <  ,  t  e  J0,  (4.3) 

where  a  is  defined  by  (3.10)  and  y  by  condition  (C). 

Proposition  3.6  provides  a  clue  how  to  proceed  because  it  suggests  that  the 
parametrization  {V,A,zQ}  of  the  original  curve  may  also  be  suitable  for  the 
approximate  curve.  With  this  in  mind,  we  define  again  an  affine  mapping 
B:  IR1  x  Y  +  X  by  B(t,y)  =  xQ  +  tzQ  +  Ay,  t  e  IR1 ,  ye  Y,  and  with  it  the 


operator 


I 


Hh:  bI'^EJcR1  x  y  -  Y,  ^(t.y)  =  Fh(B(t,y)),  (t,y)  e  B(_1  }(E). 


(4.4) 


Clearly,  Hh  is  of  class  Cr  with  derivatives 


DtHh(t,y)  =  DFh(xo*tz0*A,)zo  =  (I-P„)Qzo  *  PhDF<  VV^’V 


D  Hh(t,y)  ■  DFh(xo+tzo+Ay)A  -  ( I-Ph )QA  +  PhOF(x0+tz0*Ay)A. 


The  following  three  leumas  show  that  Hh  satisfies  the  conditions  of  Theorem 


Lemma  4.2:  For  sufficiently  small  h  and  for  each  t  e  J0,  DyHh(t,n(t) )  is 
an  isomorphism  of  Y  and  we  have  sup{ [ | DyHh(t ,n(t) )-1 | | :  t  e  JQ}  <  yQ,  where 
yQ  is  independent  of  h. 

Proof:  The  proof  proceeds  along  the  lines  of  Proposition  3.5.  In  brief,  by 

A 

condition  (C),  DyHh(0,0)  =  DFh{xQ)A  is  an  isomorphism  of  Y  with 

| | DyHh (0 .0) -1 | |  <  1/y.  From  this,  it  follows  that  DyHh(t,n(t))  is  an  iso¬ 
morphism  for  t  e  JQ  and  now  (4.3)  can  be  used  to  complete  the  proof. 


Lenina  4.3:  For  sufficiently  small  h,  we  have  sup{ | | DtH^(t,n(t) ) | | :  t  e  JQ}  <  yj 
where  y-|  is  independent  of  h. 

Proof :  This  estimate  follows  directly  from  (3.10)  and  the  continuity  of 
DF(x0  +  tz  +  An(t))  for  t  in  the  compact  interval  JQ. 


MS" 
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Laima  4.4:  lim  sup  |  |H.(t,n(t) )  |  j  =0. 

M)  teJQ  n 

Proof:  We  have 

sup  |  |H.(t,n(t))|  |  <  1 1  (I-Ph)Qx  ||  +  M|  |  (I-P  )Qz  ||  +  sup  |  |U-Ph)QAn(t)  1 1 , 
teJ  n  -no  n  0  teJ  n 

o  o 

where  M  =  max  |t|.  By  (3.3),  it  follows  that  (I-P,)Qxrt  -*■  0,  (I-P.  )Qz  0, 

teJ  n  o  n  o 

o 

and  (I-Ph)QAn(t)  -*■  0  pointwise  for  t  e  JQ.  From  (3.10)  and  the  continuity 
of  n'(t),  we  obtain,  moreover,  that  ( I-Ph)QAn(t)  -+  0  uniformly  on  JQ,  which 
implies  the  result. 

If  we  assume  now  that  OF  satisfies  a  Lipschitz  condition,  then  we  are  in 
a  position  to  apply  Theorem  4.1.  More  specifically,  suppose  that  OF  is 
lipschitz-continuous  on  bounded  subsets,  that  is,  that 

1 1 DF(x1  )-0F(x2) |  |  <  £(p)  1 1 x1  -x2 1  |  ,  x]tx2  e  E,  x?  e  Bfx^p),  (4.5) 

with  some  monotonical ly  increasing  function  £:  -*-Ir|.  Then  we  obtain  the 

following  theorem: 

Theorem  4.5:  Suppose  that  the  conditions  of  Theorem  3.7  hold  and  that  DF 
satisfies  (4.5).  By  Theorem  2.3,  the  original  problem  (2.2)  has  a  solution 
segment  x:  J  <=■  F1  -*■  E  defined  on  an  open  interval  J  c  IR1 ,  0  e  J.  Then 

there  exists  a  compact  subinterval  Jqc  j,  o  e  JQ,  such  that,  for  sufficiently 
small  h,  the  local  approximate  problem  (3.5)  has  a  solution  segment  x^:  JQ  -*•  Eh 
and 
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| |x(t)-xh{t)| |  <  C||(I-Ph)Qx(t)||,  t  e  JQi  (4.6) 

where  C  is  independent  of  h  and  t. 

Proof:  By  applying  Theorem  4.1,  we  obtain  the  existence  of  a  Cr-function 
V  +  Y  such  that  Hh(t,nh(t))  =0  for  t  e  JQ  and 

I  |nh(t)-n(t)| |  <  C0||Hh(t,n(t))|  j,  t  e  JQ,  (4.7) 

where  CQ  depends  only  on  yQ  and  and  hence  is  independent  of  h  and 
t.  With  xh(t)  =  xQ  +  tz0  +  Anh(t),  we  have  Fh(xh(t))  =  Hh(t,nhlt))  =  0  for 
t  e  JQ.  Proposition  3.1  (i)  shows  that  xh(t)  e  Eh  and  Fh(xh(t))  =  yQh; 
that  is,  x^:  JQ  ■*"  ^h  a  s°lu^^on  segment  of  the  local  approximate  problem 
(3.5).  Finally,  the  estimate  (4.6)  follows  from  (4.7)  and 

I |x(t)-xh(t) | |  =  | |A(n(t)-nh(t))| |  <  ||A||  C0||Hh(t,n(t))|| 

*  C0I|A||  ||(I-Ph)Qx(t)||,  t  e  JQ. 

It  is  worth  noting  that  Lemma  4.4  implies  a  uniformity  in  this  estimate. 

In  other  words,  we  have  the  following  result: 

Corollary  4.6:  Under  the  conditions  of  Theorem  4.5,  we  have 

x(t)  =  lim  x.  (t)  uniformly  for  t  e  J  . 

M)  h  0 


V  .. 


(4.8) 


The  remaining  three  sections  of  this  paper  are  devoted  to  examples  of 
our  general  theory.  We  begin  by  showing  that  for  mildly  nonlinear  operators 
our  assumptions,  and  especially  condition  (C),  are  naturally  satisfied.  This 
covers  the  case  considered  in  [5]. 

As  noted  at  the  beginning  of  Section  3,  in  many  applications  we  find  that 
X  =  Y  xjR1.  Let  Q  e  L( X, Y)  be  the  natural  projection 

Q(u,A)  =  u,  x  =  (u,A)  e  X.  (5.1) 

/V  A  A 

Let  X  be  another  Banach-space,  K  e  L(X,Y)  a  compact  mapping,  G:  X  ->■  X  an 
operator  of  class  Cr,  r  >_  1 ,  and  define 

F:X-*Y,  F(x)  =  Qx  +  KG(x),  x  e  X.  (5.2) 

With  x  =  (u,A),  our  problem  (2.2)  then  takes  the  familiar  form 


u  +  KG(u,A)  =  y. 


(5.3) 


Many  elliptic  boundary-value  problems  can  be  written  in  the  form  (5.3). 
For  example,  consider  the  problem 


Lu  +  g(u,A)  =  -  I  gjr-  (a,j(5)  ^7)  +  g(u,A)  =0  in  fl, 
i  ,j=l  i  J 


(5.4) 


u  =  0  on  3fi, 


where  Q  is  a  suitable  bounded  domain  in  IR  ,  the  coefficients  a.-  are 

*  J 

sufficiently  smooth,  and  the  linear  part  L  is  strongly  elliptic.  A  weak 
formulation  of  (5.4)  is 
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/ 


•  *  T 

1,J  =  1 


3ij(^  3^  3^  +  9(U’X)W 

J  ' 


d£  =  0,  we 


tyn). 


(5.5) 


Under  appropriate  regularity  conditions,  any  weak  solution  u  e  H^ft)  of 
(5.5)  turns  out  to  be  a  strong  solution  of  (5.4).  In  order  to  write  (5.5) 
in  the  form  (5.3),  suppose  that  g  defines  a  Cr-mapping,  r  >_  1 , 

G:  H^ft)  x  1R1  H-1(ft).  For  the  strongly  elliptic  bilinear  form 


a(u,w) 


1 


n 

J, 


„  /  f  \  3u  3w  jr 

ij(C)  35j  3C1 


u,w  e 


iyn). 


(5.6) 


let  K:  H’^fft)  -*•  Hg(ft)  be  the  linear  operator  defined  by 

a(Ku,w)  =  (u.w )Q,  u  e  H-1(ft),  w  e  H^(ft),  (5.7) 

2  2  1 

where  (u,w)Q  is  the  inner  product  on  L  (ft).  Then  K:  L  (ft)  -*■  HQ(ft)  is 
compact  and  the  weak  formulation  (5.5)  is  equivalent  to  (5.3)  for  yQ  =  0. 

As  in  Section  3,  we  introduce  a  family  {P^}  of  projections  Ph  e  L(Y) 
of  finite  rank  such  that  (3.3)  holds  and  define  the  subspaces  Y^  =  P^Y, 

Xh  =  Yh  x]r1,  and  Zh  =  Yh  *  {0>’  h  > 

A 

With  as  defined  in  (3.7),  the  following  lemma  is  the  key  to  the 
applicability  of  our  results  to  (5.3): 


Lemma  5.1 :  For  the  problem  (5.3)  and  for  any  fixed  xQ  =  (u0,AQ)  e  E,  we  have 

lim  DF.(x  )  =  DF(xn)  in  the  uniform  operator  topology  on  L(X,Y). 
h-»0  n  0  0 

A 

Proof:  Evidently,  the  definitions  of  the  operators  F  and  F^  imply  that 
DF(xq)  -  0F^(xo)  =  (I-Ph)K0G(xQ) ,  and  the  compactness  of  K  ensures  that 
(I-Ph)KDG(x0)  ■*  0  in  the  uniform  operator  topology  on  L(X,Y). 


Next,  let  {V,A,zq}  be  a  local  parametrization  of  M(yQ)  at  a  point 
xo  =  <V*o>  e  Then  we  can  ver^y  the  following  key  fact: 

Lemma  5.2:  For  problem  (5.3),  condition  (C)  is  satisfied. 

Proof:  Once  again  we  proceed  as  in  the  proof  of  Proposition  3.5.  Clearly, 

DF(xq)A  is  an  isomorphism  of  Y,  and  hence  the  result  is  a  direct  consequence 

of  the  fact  that,  by  Lemma  5.1,  lim  DF.  (x)A  =  DF(x  )A  in  the  uniform  operator 

h-0  h  0  0 

topology  on  L(Y) . 

With  this,  we  have  shown  that,  for  mildly  nonlinear  operators,  our  basic 
assumptions  indeed  are  satisfied  and  hence  that  Theorems  3.7  and  4.5  apply.  We 
note,  in  particular,  the  form  of  the  estimate  (4.6)  in  this  case.  If,  as  in 
Theorem  4.5,  x(t)  =  (u(t),A(t)),  t  e  JQ,  and  xh(t)  =  ( uh ( t ) ,Ah ( t) ) ,  t  e  JQ, 
denote  the  points  on  the  solution  segments  of  the  original  and  approximate 
problems,  respectively,  then  (4.6)  becomes 

||u(t)-uh(t)||  +  |A(t)-\h(t)|  <  C||(I-Ph)u(t)j|,  t  e  J0.  (5.8) 

All  points  on  regular  solution  curves  may  be  classified  as  being  either 
nonsingular  points  or  limit  points.  It  is  interesting  to  see  what  types  of 
local  parametrizations  correspond  to  these  two  types  of  points  as  considered 
in  the  first  two  parts  of  [5]. 

In  defining  nonsingular  points  and  limit  points,  there  is  no  need  to 
restrict  ourselves  to  mildly  nonlinear  mappings. 

Definition  5.3:  Let  (A)  be  given  and  suppose  that  a  splitting  (3.1)  is  available 


together  with  the  operator  Q  of  (3.2).  Then  —  with  a  mild  abuse  of 
notation  --  a  point  xQ  =  (u0,XQ)  e  E,  uQ  e  Z,  XQ  e  T,  is  nonsingular  if 
DuF(xo)  is  an  isomorphism  of  Z  onto  Y.  The  point  is  a  1 imit  point  (with 
respect  to  the  direction  of  T)  if  DuF(xQ)  has  a  one-dimensional  null  space 
in  Z  and 

DxF(xQ)  4  rge  DuF(xq).  (5.9) 

The  regularity  of  such  points  is  the  content  of  the  following  result: 

Proposition  5.4:  Let  the  information  of  Definition  5.3  be  given.  If  xq  £  E 
is  either  a  nonsingular  point  or  a  limit  point,  where  DF(xq)  e  L(X,Y)  is  a 
Fredholm  operator  of  index  1,  then  xQ  e  R(F). 

Proof :  If  xQ  is  a  nonsingular  point,  then,  since  DuF( xQ )  maps  onto  Y, 
we  have  rge  DF(xq)  =  Y,  and,  since  DF(xQ)  is  a  Fredholm  mapping  of  index  1, 
it  follows  that  dim  ker  DF(xQ)  =  1  and  hence  that  *Q  e  R(F).  If  xq  is  a 
limit  point,  then  the  formula  DF(xJ(w,y)  =  D  F(xjw  +  yD.F(x  )  for  (w,y)  e  X, 
together  with  (5.9),  implies  that  ker  DF(xQ )  =  ker  DUF( xQ )  *  {01,  and  thus 
dim  ker  DF(xQ)  =  dim  ker  duF(xq)  =  1.  Since  DF(xQ)  is  a  Fredholm  operator 
of  index  1,  it  follows  that  rge  DF(xQ)  =  Y  and  again  that  xQ  e  R(F). 

We  return  again  to  the  mildly  nonlinear  case  (5.2).  Clearly,  since  Q  is 
a  Fredholm  mapping  of  index  1,  so  is  DF(x)  =  Q  +  KDG(x)  for  any  x  e  X.  Thus, 
Proposition  5.4  applies.  If  xQ  is  a  nonsingular  point,  then,  because  of 
(Y  x  {0})nker  DF(xQ)  =  (0),  we  may  choose  as  a  local  parametrization  the  subspace 
V  =  Y  x  {0},  the  isomorphism  A:  Y  -*•  V  defined  by  Aw  =  (w,0),  w  e  Y,  and 
zQ  =  (0.1).  This  parametrization  amounts  to  the  selection  of  X  as  the  parameter 
and  corresponds  in  essence  to  the  choice  in  the  first  part  of  [5].  The  estimate 


(5.8)  now  reduces  to 


||u(X)-uhU)||  <  C||(I-Ph)u(X)||,  X  e  JQ.  (5.10) 

If  e  X  is  a  limit  poirt,  then  with  u.  e  1  such  that  D  F(x  )u  =  0, 
o  r  o  u  o  o 

llu0ll  =  1 »  and  Y0  =  span  {u0},  Y-j  =  rge  DuF(xQ),  it  follows  tnat  Y  =  Y-j  ©  Y  , 
OuF(xo)  is  an  isomorphism  of  Y-j ,  and  ker  DF(xQ)  =  YQ  *  {0}.  Since 
(Y1  xJI1)  A  ker  DF(xQ)  =  {0},  we  may  use  the  local  parametrizat'on  specified 
by  V  =  Yj  x  the  isomorphism  A:  Y  -*•  V  defined  by  Aw  =  (w-j  ,to)  for 
w  =  wi  +  tQuo  e  Y,  w-j  e  Y-j ,  and  zo  =  (uQ,0).  This  corresponds  exactly  to 
the  approach  in  the  second  part  of  [5]. 


6.  More  General  Nonlinear  Problems 

The  techniques  and  results  presented  here  provide  a  unification  which  is 
certainly  of  interest  in  itself,  but  their  real  value  should  derive  from  the 
fact  that  the  results  are  applicable  to  a  wider  class  of  problems  than  the 
mildly  nonlinear  problems  to  which  other  known  approaches  appear  to  be  restricted. 
The  full  extent  of  this  applicability  is  still  under  study,  but  our  results 
certainly  cover  problems  which  to  our  knowledge  could  not  be  handled  by  other 
methods.  A  model  problem  of  this  type  is  discussed  in  this  section. 

The  flexibility  in  our  formulation  lies  in  the  freedom  of  choice  of  the 
operator  Q  of  (3.2).  In  the  case  of  the  mildly  nonlinear  operators  discussed 
in  the  previous  section,  Q  may  be  chosen  as  the  natural  projection  (5.1).  For 
more  general  operators,  Q  has  to  reflect  more  specifically  some  of  the 
properties  of  the  problem.  The  flexibility  in  the  choice  of  Q  may  even  be 
important  in  the  case  of  mildly  nonlinear  boundary- value  problems  such  as  (5.4). 


Although  (5.4)  can  be  transformed  into  a  mildly  nonlinear  problem,  the 
techniques  to  be  discussed  below  also  allow  for  a  direct  use  of  the  strong 
form  (5.4). 

As  our  model  problem,  we  consider  here 


arc  tan  (^|)  +  g(u,A)  =0,  0  <  £  <  1, 
u(0)  =  u(l)  =  0, 


(6.1) 
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where  g:  IR  1R  is  a  sufficiently  smooth  function.  A  problem  of  this  type 
was  used  in  [3J  as  an  example-case  for  numerical  experimentation  with  problems 
concerning  the  deformation  of  nonlinear  rods. 

An  attempt  to  analyze  (6.1)  by  means  of  a  weak  formulation  as  in  the 
previous  section  will  run  into  difficulties  because  the  resulting  operator  F 
does  not  satisfy  the  required  continuity  and  Lipschitz  conditions.  This  is  due 
to  the  nonlinearity  arising  from  the  arctan  function.  In  order  to  work  with  (6.1) 
directly,  let  Z  =  C^[0,1]  and  Y  =  C°[0,1],  and  consider  the  operator 


F:  X  =  Z  *  IR  -*■  Y,  F(x)  *  a [arctan(du)]  +  g(u,A), 


where,  for  abbreviation,  we  use  du  =  u.  The  derivative  of 


x  =  (u,A)  e  X, 

(6.2) 

F  is  given  by 


DF(x)(w,y)  =  d 


1 


1  +  (du)‘ 


dw 


+  Dug(u,A)w  +  pDxg(u,A),  x  =  (u,A)  e  X, 


(w,y)  e  X.  (6.3) 


It  is  a  straightforward  (but  somewhat  tedious)  matter  to  show  that  F  is  in- 


finitely  differentiable  and  DF  is  Lipschitz-continuous  on  bounded  subsets 
in  the  sense  of  (4.5). 

For  G:  X  -*•  Y,  G(x)  =  g(u,A),  x  =  (u,X)  e  X,  it  follows  from  Ascoli's 
Theorem  that,  for  fixed  xQ  =  (uQ ,Aq)  e  X,  DG(xQ)  e  L(X,Y)  is  compact.  In 
view  of  this  and  (6.3),  we  now  define 

(6.4) 


(6.5) 


q  e  L(X,Y),  Q(w,p)  -  d[" - - - 5-  dw  |,  (w,p)  X. 

1’  +  (duo>  J 


It  is  easy  to  prove  the  following  two  properties  of  Q: 


(i)  ker  Q  =  {0}  xR1, 


(ii)  Q|Z  x  {0}  is  an  isomorphism  from  Z  *  {0}  onto  Y. 


In  terms  of  Q  and  DG  the  derivative  (6.3)  of  F  at  xQ  has  the  form 

0F(xo)  =  Q  +  DG(xq).  (6.6) 

Now  proceed  as  before  and  let  (Ph)  be  a  family  of  projections  Ph  e  L(Y) 
of  finite  rank  and  define  the  subspaces  (3.4).  With  Fh  as  in  (3.7),  we  have 
DF(xq)  -  DFh(xQ)  =  (I-Ph)DG(xQ),  and  as  in  Section  5  we  can  show  that 

A 

lim  DF.  (x. )  =  DF(x„)  in  the  uniform  operator  topology  on  L(X,Y).  It  follows 
M)  "  0  0 

in  turn  that,  for  any  local  parametrization  {V,A,zQ}  of  M(0)  at  a  point 
xo  =  ^uo’*o^  e  the  condition  (C)  holds. 

Clearly  now,  for  the  mapping  F  of  (6.2)  and  with  Q  defined  by  (6.4), 
our  results  apply  to  the  problem  (6.1).  If  x(t)  =  (u(t),A(t)),  t  e  0Q,  and 
xh(t)  =  (uh(t),Xj)(t)),  t  e  JQ,  once  again  denote  the  points  on  the  solution 
segments  of  the  original  problem  and  the  approximate  problem,  respectively,  then 
the  estimate  (4.6)  takes  here  the  form 
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|u(t)-uh(t)||  +  |A(t)-Ah(t)|  <  C 

C2 


(I-Ph)  d 


1  +  (duQ)‘ 


du(tx 


*  t£V 


(6.7) 


The  mapping  F  of  (6.2)  serves  only  as  an  example  of  the  class  of  non¬ 
linear  operators  to  which  the  approach  in  this  section  may  be  applied.  In 
principle,  let  F  be  any  mapping  of  the  form 

F:  X  =  Z  x  IR1  +  Y,  F{ x )  =  N(x)  +  G(x),  x  e  X.  (6.8) 

If  xQ  e  X  is  a  point  for  which  Q  =  DN(xQ)  has  the  properties  (6.5)  and 
DG(xQ)  is  compact,  then  clearly  the  results  apply. 

We  note  in  passing  that  the  general  discussion  of  nonsingular  points  and 
limit  points  at  the  end  of  Section  5  and  the  corresponding  parametrizations 
carries  over  to  mappings  of  the  form  (6.8). 


7.  The  Reduced-Basis  Technique 

Our  final  example  concerns  the  so-called  reduced-basis  technique  which  has 
been  receiving  increased  attention  in  the  engineering  literature  (e.g.,  see  the 
survey  [11]  where  other  references  are  also  given).  The  technique  arises  in  the 
context  of  a  standard  continuation  procedure  for  nonlinear  problems  and  is  usually 
viewed  as  a  method  for  effecting  considerable  reductions  in  the  size  of  the 
systems  of  equations  obtained  from  finite-element  approximations.  Alternately,  as 
discussed  in  [6],  the  technique  corresponds  to  the  construction  of  a  set  of  finite- 
element  basis  functions  which  closely  reflects  the  properties  of  the  solution 
segment  that  is  to  be  approximated.  In  this  form,  the  reduced-basis  technique  turns 
out  to  fit  very  naturally  into  the  general  setting  of  this  paper. 
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As  noted  earlier,  condition  (C)  establishes  a  certain  compatibility 
between  the  discretization  and  local  parametrization.  In  the  previous  two 
sections,  we  saw  that  its  validity  can  be  ensured  for  certain  classes  of 
operators  without  restricting  the  choice  of  the  discretization  and  parametrization. 
Conversely,  it  is  possible  to  enforce  condition  (C)  without  placing  any  further 
conditions  on  the  mapping  F  by  selecting  certain  types  of  discretizations  and 
parametrizations  inherent  to  the  problem.  The  reduced-basis  technique  represents 
an  example  of  this  latter  approach. 

As  usual,  suppose  that  the  information  (A)  is  given  and,  for  ease  of 
notation,  assume  that  the  mapping  F  is  defined  on  all  of  X.  For  a  fixed 
yQ  e  F(R(F) )  and  xQ  e  M(yQ )  we  consider  the  mapping 

Fo:  x  -  Y.  Fq(x)  =  F(x0+x),  x  e  X,  (7.1) 

and  note  that  0  e  M0(yQ)  =  Fg_1^(yo)  A  R(FQ).  Now  our  procedure  is  applied 

to  F„  rather  than  F. 

o 

For  a  splitting  (3.1)  of  X,  we  choose  X  =  Z  ©  T,  where  T  =  ker  DF(xq) 
and  Z  is  any  complementary  subspace.  Then  Q  =  DF(xq)  e  L(X,Y)  satisfies 
ker  Q  =  T,  and  Q|Z  maps  Z  isomorphical ly  onto  Y.  As  usual,  let  a  family 
{P^}  of  finite-rank  projections  of  Y  be  given  for  which  (3.3)  holds  and 
define  the  subspaces  (3.4).  For  a  local  parametrization  {V,A,zQ},  we  now 
choose  V  =  Z,  A  =  Ay  =  (Q|Z)’\  and  any  nonzero  vector  zQ  t  T.  Then  we  have 

DF0h(°)A  =  (I-Ph)QA  +  PhDFQ(0)A 

-  (I-Ph)DF(x0)A  +  PhDF(xo)A 


•  DF(x0)Ay  =  I, 


where  I  is  the  identity  on  Y.  In  other  words,  condition  (C)  holds  trivially, 
regardless  of  the  nature  of  F. 

Now  Theorem  4.5  applied  to  Fq  ensures  the  existence  of  an  approximate 
solution  segment  xQh(t)  =  0  +  tzQ  +  Anh(t)  e  Xh,  t  e  J  ,  satisfying 

PhFo(xoh(t))  =  yoh*  teV  before,  x*(t)  =  xQ  +  xQh(t)  lies  in  xQ  +  Xh 

and  satisfies 

fhF<xhR(t»  ■  ^oh-  teV  l7'2» 

From  the  estimate  (4.6),  we  see  that  xQh(0 )  =  0  and  thus  xP(0)  =  xQ.  By 

differentiation  of  (7.2)  with  respect  to  t  at  t  =  0,  it  follows  that 

0  =  PhDF(xo)(zo  +  An£(0))  =  Phn^(0).  (7.3) 

Since  zq  and  tzQ  +  Anh(t)  are  in  Yh>  the  same  holds  for  nh(t),  t  e  JQ, 

whence  (7.3)  implies  that  n^(0)  =  ^  or  xP(0)  =  zQ* 

D 

Hence,  the  solution  segment  xh(t),  t  e  JQ,  passes  through  xQ  and  has 
the  same  tangent  direction  as  the  original  curve.  By  appropriate  choice  of  the 
space  Z  and  the  projections  P^,  it  is  easy  to  ensure  that  higher  derivatives 
of  the  original  curve  and  the  approximate  curve  also  agree  at  xQ.  (Note, 

however,  that  the  validity  of  condition  (C)  requires  only  that  ker  DF(xq)  c 

and  not  any  additional  matching  of  higher  derivatives  at  x  . )  This  is  precisely 
the  approach  of  the  reduced-basis  technique;  that  is,  the  method  uses  the  first 
few  directions  of  the  moving  frame  of  the  original  solution  curve  WQ(y0)  at 
the  origin  as  the  basis  vectors  for  the  approximating  subspaces  X^.  For  some 
details  and  an  error  assessment, we  refer  to  [6].  As  the  mentioned  survey  [11] 
already  indicates,  the  method  is  certainly  beginning  to  prove  itself  very 
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effective  for  the  solution  of  various  geometrically  nonlinear  structural 
probl ems . 
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